1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 4000000 iteration;
- 3000000 burn-in;
- 1000 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1 1.01
## LDevsum 1 1.02
## dh0 1 1.00
## dl0 1 1.00
## dl1 1 1.00
## dl2 1 1.00
##
## Multivariate psrf
##
## 1.01
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1364.55658 | 22733.4424 |
| DIC3 | 1246.39509 | 22682.1610 |
| PWAIC | 83.75423 | 292.0923 |
| WAIC | 1307.00071 | 22709.9272 |
2. Separate Modeling of High-Level
2.1 Trace plots for convergence check
The current MCMC setting is:
- 4000000 iteration;
- 3000000 burn-in;
- 1000 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 0.999 0.999
## dh0 0.999 0.999
##
## Multivariate psrf
##
## 0.999
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| H0b | |
|---|---|
| DIC | 1392.95198 |
| DIC3 | 1260.28425 |
| PWAIC | 94.66716 |
| WAIC | 1332.37089 |
3. Separate Modeling for Low-level
3.1 Trace plots for convergence check
The current MCMC setting is:
- 4000000 iteration;
- 3000000 burn-in;
- 1000 thinning.
3.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## LDevsum 0.999 1.00
## dl0 1.001 1.01
## dl1 1.004 1.02
## dl2 0.999 1.00
##
## Multivariate psrf
##
## 1
3.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
3.4 WAIC results
| L2 | |
|---|---|
| DIC | 22588.9643 |
| DIC3 | 22659.4118 |
| PWAIC | 322.1958 |
| WAIC | 22695.5826 |